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Abstract 



We present variational calculations of the one-body density matrices and 
momentum distributions for ^He-'^He mixtures in the zero temperature limit, 
in the framework of the correlated basis functions theory. The ground-state 
wave function contains two- and three-body correlations and the matrix ele- 
ments are computed by (Fermi)Hypernetted Chain techniques. The depen- 
dence on the ^He concentration (3:3) of the ^He condensate fraction (^0^^) 
and of the ^He pole strength {Zp) is studied along the P = isobar. At low 
^He concentration, the computed ^He condensate fraction is not significantly 
affected by the ^He statistics. Despite of the low X3 values, Zp is found to be 
quite smaller than that of the corresponding pure ^He because of the strong 
^He-^He correlations and of the overall, large total density p. A small increase 
of riQ^* along X3 is found, which is mainly due to the decrease of p respect to 
the pure "^He phase. 
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The momentum distributions (MD) of atoms in quantum hquids is a challenging problem 
of fundamental interest .011 They provide essential information on the correlations present in 
the system, which do not show up explicitly in other quantities. In the past years, accurate 
inelastic neutron scattering experiments have allowed for studying several aspects of the 
momentum distribution in helium liquids, ^He,i'i ^Hei and ^He-^He mixtures.i'0 However, 
a clean extraction of information on the Helium MD's is somehow tampered by the need of 
a sound theoretical understanding of the final state effects in the analysis of the dynamic 
structure function, even at high momentum transfers. 

The theoretical methods to evaluate momentum distributions of many-body interact- 
ing, dense systems at zero temperature have also made a significant progress in recent 
years.lll At present, there are results for the pure Helium phases obtained within different 
many-body techniques, i.e., variational theory (using either integral equations!'! or Monte 
Carlo methods0) and almost exact stochastic methods as Green's Function Monte Carlo 
(GFMC)0li or Path-Integral Monte Carlo (PIMC ).i 

The MD's of liquid ^He (^He) are influenced by the Bose (Fermi) statistics of the atoms. 
The macroscopic occupation of the zero momentum state, as given by the condensate fraction 
n^^ , characterizes the momentum distribution of bosonic, liquid ^He and it is strictly linked 
to its superfiuid behavior. On the other hand, the discontinuity Zp at the Fermi momentum 
kp is a. characteristic of the ^He system when it is studied as a normal Fermi liquid. 

In this paper we consider the interesting case of isotopic ^He-^He mixtures where, due to 
its fermion -boson nature, both quantities Zp and n^^ are simultaneously present. Recent 
neutron scattering experiments on Helium mixtures at high momentum transfer^'! give 
additional motivations to undertake a microscopic, theoretical study of their momentum 
distributions and one-body density matrices. Special emphasis will be devoted to the de- 
pendence on the ^He concentration, X3, of the single-particle kinetic energies of the isotopes 
and of Zpand 

The investigation is carried on in the framework of the variational approach. The 
trial wave function for the mixture contains two-body (Jastrow) and triplet correlations. 



This type of correlated wave function has been useful in effectively studying the pure 
phases.™ Two of usill (A. P. and A.F.) derived the hypernetted and Fermi hypernet- 
ted chain (HNC/FHNC) equations for the momentum distributions of the mixtures using 
trial wave functions with only pair correlations. Numerical applications were carried out 
in the HNC/FHNC/0 approximation, i.e., neglecting the elementary diagrams. A prelimi- 
nary study of the elementary diagrams for a Jastrow trial wave function was performed^ 
by generalizing the scaling approximation proposed for pure phases.i'i Also available are 
variational Monte Carlo (VMC) calculations0 with similar correlations of the analytical 
McMillan type. The studies of the mixture has been recently complemented with varia- 
tional calculations concerning the energy and stability of the ground state, with path 
integral Monte Carlo (PIMC) analysi£ll and with microscopic correlated basis functions 
(CBF) estimates of the inelastic neutron scattering cross sections both at intermediateil 
and hig hi momentum transfers. 

The paper is organized as follows: in the second section we will shortly present the 
HNC/FHNC theory to calculate n{k) for mixtures described by correlated wave functions 
containing two- and three-body correlations. The treatment of the elementary diagrams in 
the so called scaling approximation is discussed in some details in the second part of the 
section. Results for n^'^\k), n^^\k) and for the one-body density matrices are presented in 
Section II, together with a critical discussion of the discrepancies with the available analysis 
of the deep inelastic neutron scattering measurements on mixtures, which (in contrast with 
our results) point to a large enhancement of the ^He condensate fraction. 

I. HNC/FHNC EQUATIONS FOR THE MOMENTUM DISTRIBUTION OF 

^'HE-'^HE MIXTURES 

The one-body density matrices p^^'^ri, r[) (a = 3, 4) for a homogeneous, isotopic mixture 
of N3 ^He atoms and A''4 ^He atoms, described by a ground-state wave function "^{1, . . . , N4 + 
N3) are defined as 
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iV, / ...,N, + Ns)^iV^, ■ ■ ■ , iV4 + iV3)c;r2 . . . dr^,+N, 

^ ^""''""'^ Pa /|vl/(l,...,iV4 + iV3)|2ciri...rfr^,+jV3 • ^' 

In homogeneous mixtures, with constant particle densities pa = Na/N, p'^"^(ri, r'^) = 
p(")(r), with r =| ri — r'^ |. p*^"^(r)'s satisfy the normahzation conditions Uap^^'H^) — 1; 
z/q, being the spin degeneracy (u^ = 1, = 2). Notice that in the definition of p^^\r) the 
spin variables have not been explicitly written. We will henceforth omit the subindex in the 
degeneracy factor and assume that it always refers to '^He. 

The momentum distribution of the a component, or rather the occupation probability 
for single-particle states with momentum k and given spin projection, can be obtained as 
the Fourier transform of the corresponding density matrix, 

n(°) (A;) = 5a4P4ni,'^ (27r)=^5(k) + p„ j rfr exp(zk • r) [p(") (r) - 6^,n^^^] , (2) 

where n^^^ = p'^^^oo) is the ^He condensate fraction, i.e., the fraction of ^He particles in the 
zero momentum state. 

The ground state of the mixture is well described by a generalization of the correlated 
wave function used in the pure phases: 

vi>(i,...,iV4+iV3)= n n /^°'^H^«,j/3) n f^'''^''\^a,J(^,k,)<p{l,...,N,). (3) 

a</3<7=3,4 ia<j/3 ia<jf3<kj 

0(1, . . . , N^) is the Slater determinant of plane waves corresponding to the Fermi component 
of the mixture, and f^°''^\iai jp) ( f^"''^'"'K'i'a,jf3, ky)) are the 2 (3)-body correlation functions 
involving 2 (3) particles of types a, (3 (a, /3, 7), respectively. Similar trial wave functions have 
been used in previous works to study the structure and energetic ground-state properties of 
=^He-^He mixtures. 

A cluster analysis of p(")(r) in powers of w^"'^) = /("'^^ - 1, /i^^'^^) = [f^'^'^'^f - 1, 
^(«,/3,7) = /(«,/3,7) _ 1 and = [/("'/5,7)]2 _ 1^ as that carried out in the pure phases,!!!! 

gives the following structural decomposition for 

p(")(r) =nS,")Ar(°)(r) , (4) 
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where massive re-summations of the diagrams, as defined in Refs. p,^,|T^,^, may be per- 
formed in practice by using HNC/FHNC techniques.§@'il 
The strength factor ng"^ is given by 



(a) 



exp[2r(r)-ri")] (5) 

and 

iV-(r) = [5,4 + S^silKkpr) - iV^Kr) - i?iKr))] exp [iV;("J(r) + E("J(r)] (6) 

sums up all the irreducible diagrams with external points la and 1^. In Eq.(|D, /(x) = 
3ji{x)/x is the Slater function and kp = (Gn'^p/i'Y^^ is the ^He Fermi momentum. 

The functions Nl^\r) and E'^{r) are the sums of the nodal and elementary diagrams 
contributions, respectively. The evaluation of the nodal functions N^y\r)^ in the context 
of the HNC/FHNC approach, is discussed in the Appendix, also containing the explicit 
expressions of the F^"^ factors. 

The momentum distributions are computed via the density matrices by Eq.(0). We thus 

get 

n^'\k) = (27r)V4n^)5(k) + p,ni^^ J rfrexp[zk ■ r](exp[iV(t)(r) + Ej^^ir)] - 1) , (7) 



and 



where 



n(3)(A;) = nP[n,{k) + Q^kp - k)na{k)] , (8) 



= 1 - X,, + 2X^^, + (9) 



and 



nc{k) = -—^ - P3 / cir exp[^k • r] {(exp[iV£)(r) + E^l{r)] - 1) (10) 

t — ^cc 

x(-/(A;^r)/z/ + A^(l(r) + E(fiJr)) + E(l(r)} . 



Xyc = Qyc — Nyc + for y = Uc, c and Xxy{k) stands for the Fourier transform 

X^y{k) = p^ j dre'^'-'X^yir) (11) 

The strength factor n^^ is the asymptotic value of the ^He one-body density matrix, 
pi^)(r 00) = nl^^ and corresponds to the ^He condensate fraction. The decomposition 
of rS^\k) in a continuous (nc(/i;)) and a discontinuous {nd{k)) piece exphcitly hnks the 
discontinuity of n^^\k) at kp, Zp, to Udikp) by 

Zp = nfna{kp) . (12) 



A. Scaling approximation for the elementary diagrams 

The HNC/FHNC equations can be solved once a given prescription for the contributions 
of the elementary diagrams has been given. However, as no exact method to compute them 
is presently known, at least in the frame of the integral equations, one has to resort to some 
approximation. Among the available schemeJ^Zhil we have chosen the scaling approximation 
(SA), developed for both the energy and the one-body density matrix of pure phases,ili01li 
and satisfactorily reproducing VMC calculations. Although the number of elementary di- 
agrams in the mixture is much larger, it is straightforward to generalize the pure phases 
scaling approximation to our case. 

The SA is based on the evaluation of the 4-points elementary diagrams constructed with 
the combinations of the distribution functions gi^'^\r) allowed by the diagrammatic rules 
and it has already been used in the calculation of the energy and of the static structure 
functions of the mixture.0 The elementary diagrams are approximated by 

E£'^)(r) = E{r) ,Ei'^/\r) = 0, G {3,4}, xy = [de,ee,cc] , (13) 



where 



E{r) = (1 + s)EW(r) + El^\r) . (14) 
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Ef\r) and (r) are the four-points elementary diagrams without and with explicit three- 
body correlations into their basic structure, respectively. These diagrams are constructed 
by using as internal links an averaged dressed correlation g{r) — 1, 

g{r) = xlg^'^'\r) + 2x,x,g^^'^\r) + xlg^^^'\r) , (15) 

with Xa = Pa/ P- The introduction of g{r) makes feasible the calculation of E{r) because it 
reduces drastically the high number of elementary diagrams originated by all the possible 
bonds between ^He and ^He particles. Actually, for the underlying boson-boson mixture (i.e. 
$(1, . . . , A/'g) = 1 in Eq. (D) and taking the same correlation functions between all types of 
isotopes ( average correlation approximation (AC A) ), g{r) provides the exact E^^j-{r). This 
property and the small ^He concentration in the physical region of interest (0:3 < 0.10) justify 
the use of g{r). The scaling parameter s (|l^ is determined by imposing the consistency 
between the Pandharipande-Bethe and the Jackson-Feenberg forms of the kinetic energy for 
the boson-boson mixture without triplet correlations, s is calculated for each total density 
and it is kept fixed when X3 changes. This assumption is plausible because, at low ^He 
concentrations, the statistical effects in g{r) are negligible. 

The additional elementary diagrams needed for the one-body density matrices are simi- 
larly evaluated: 

ElTir) = EUr), 4?'^^ = iyz = u;e,u;,c) (16) 

with 

EUr) = (1 + s^,)ESl{r) + 41,(r), (17) 

and 

4"J(r) = (1 + sI-^)e!^1^ + Ej^Ur), (18) 
£^SL(r) = (1 + s^^^ML, + E}^ljr). (19) 
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The average distribution function 

9M) = xl gl'fir) + 2x3X4 {g^jf\r) + g^JfHr)) + {gSf\r) + g^Jf\r)) (20) 

has been used to compute the above four-points elementary diagrams. 

Finally, the set of single external point elementary diagrams, appearing in the strength 
factors no"'' expressions, are approximated, as in the pure phases,!'! by 

i^. = (l + ^sxrf)41 + 4!l, x = cu,d. (21) 

As far as the factors related to the momentum distributions are concerned, we have 
chosen s^^d by imposing Tmd = Tjp, where Tmd is the total kinetic energy obtained by 
integrating the momentum distribution, 

Tmd = I e n^^\k) + I d\. ^ n^'^k), (22) 

and Tjp is the ground-state expectation value of the kinetic energy operator computed by 
the Jackson-Feenberg identity. Moreover, the fulfillment of the normalization conditions of 
the momentum distributions, i.e.. 



(27r)3p„ 

equivalent to z/Q,p(°^(0) = 1, requires 



dkno,{k) = 1 , (23) 



(a) 

Tig exp 



A^tnO) + i5tnO)l = 1 (24) 



A^l(0) + i?SL(0)=0. (25) 

These conditions are used to determine the remaining scaling parameters {s'^], s^^^J. 

As a matter of fact, the use for the triplet correlated wave function of the same sj^^ 
and s^^^^ parameters, as determined in the Jastrow case, produces significant deviations of 
the above normalizations from their exact values. For this reason and to ensure the correct 
normalizations of the density matrices, we have recalculated the scaling factors s^jd, s^^}^, s^^2 
and s^^^^ when the wave function contains three-body correlations, as in Ref. 
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II. RESULTS 



In this section we report results for the momentum distributions of ^He-^He liquid mix- 
tures using the Aziz potential (HFDHE2)§ for the variational determination of the ground- 
state correlations. This interaction effectively describes the equation of state of the pure 
phases-lliill The interatomic potential in isotopic mixtures is the same between any pair of 
particles. Based on this fact, we have used the average correlation approximation, ACA. 
The ACA approach, which has been carefully analyzed for the impurity problem]^ has also 
been used in the past to study finite concentration Helium mixtures.0@il The potential is 
strongly repulsive at short distances, so the correlation functions are expected to show the 
same short-range behaviors. Small differences can arise however at intermediate and large 
distances, where the interaction is weaker, because of the different masses and statistics of the 
isotopes. Nevertheless, ACA may well serve to the purpose of studying the Xs-dependence of 
the momentum distributions in the mixture. In fact, for Jastrow correlated wave functions 
we have released the ACA, allowing for different correlations in different isotopic pairs, and 
these extra variational degrees of freedom have not significantly changed our results. 

The two-body correlation function /(r) has been taken to have an analytical form, of the 
McMillan type at short distance and with enough flexibility to adjust to the optimal pure 
^He correlation behavior in the intermediate and long ranges. 



f{r) = exp 



1 fbV' 



A + Bexp 



(26) 



2\rJ 

The long range, behavior ensures the proper linear dependence of the ^He structure 
function at k ^ 0. 

The /(r) parameters at the ^He energy variational minimum, at equilibrium density 
Po = 0.365 a-^ {a = 2.556 A), are b = 1.18 a, A = 0.85, B = 1 - A, D = 3.8 A and 
r = 0.043 A~^. B and r are related to the experimental pure "^He sound velocity c and to 
the \ow-k behavior of its static structure function by 

B m^c 



T 27r2^po 



(27) 



The three-body correlation function f{rij,rik,rjk) has the parameterized form:i'illllli 



f{rij,rik,rjk) = exp 

where 



^ /=0,1 eye 



(28) 



ii{r) = {r - 5/on«)exp 



(29) 



The values of the triplet functions parameters have been taken from Ref. |lj omitting the 
small / = 2 component. 

The calculations presented here are performed at the experimental values of the density 
along the P = isobar. In this regime, the density decreases from p = pQ [x^ = 0) to 
p = 0.3582 a^^ at = 0.066, corresponding to the ^He maximum solubility. The partial 
■^He density increases from zero up to p^ = 0.0236 cr~^ in the same X3 range. So, we have 
neglected the density dependence of the variational parameters of the correlations because of 
the small variations both of the total and partial densities in the region of physical interest. 

Before presenting the results for the Helium mixtures, it is worthwhile to study the 
accuracy of the scaling approximation in the case of pure ^He. We have considered a 
correlated wave function containing McMillan two-body correlations {A = 1., B = 0., 
and b = 1.20 a in Eg . (|26|) ) and a three-body factor given by Eq. (|28|) . At po we obtain 
nQ^\jTi) = 0.078 and nQ^\jT()i) = 0.081, where the JTi (JTqi) results include triplet 
correlations contributions without (with) the / = component. The corresponding en- 
ergies are E/N{JTi) = -6.55 K and E/iV(JToi) = -6.58 K. A VMC study by one of 
the authors (J.B.), with the same trial wave functions, gives nQ'^'*(JTi)(VMC) = 0.078, 
ni'^^(JToi)(VMC) = 0.082, E/A^(JTi)(VMC) = -6.617 K and E/A^( JTqi) (VMC) = -6.625 
K. These results have been confirmed by an independent VMC calculation of Moroni,ii who 
gets n^Q^ = 0.077 and E/N = -6.604 K for the (JTi) case. 

The agreement between HNC and VMC results gives confidence in the scaling approxi- 
mation to the elementary diagrams as described in the previous section, prescribing a recal- 
culation of the scaling parameters directly associated to the momentum distribution after 
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the inclusion of the three-body correlations. Actually, if the scaling parameters in the JT 
cases are the ones determined at the Jastrow level (as in Refs. gig), we get 4'^^(JTi) = 0.064 
with a violation of the normalization conditions of ~ 15%. In addition, the / = component 
of the triplet correlation has been found to have a very small effect on both the energy and 
condensate fraction. This finding also has been confirmed by the Moroni calculations^ and 
is in contrast with that of Refs. ^|36|, where the relative change in tiq was about 25 %. Due 
to the small effect of the / = triplet correlation, we have omitted its contribution in all the 
results presented for the mixture. 

The use of the semi-optimized two-body correlation factor of Eq. (p6D and of the / = 1 
triplet correlation lowers the energy to —6.62 K and provides n'^^ = 0.082. The Euler Monte 
Carlo (EMC) result of Ref. |3^, using fully optimized two- and three-body correlations in a 
VMC scheme, is n[,'^^(EMC) = 0.087. On the other hand, the DMC results of Refs. 0,|l| are 
nf,'^^(DMC) = 0.072 and nS,'^^(DMC) = 0.084, respectively. The difference between the two 
DMC results is due to the use of an extrapolated estimator which is sensitive to the overlap 
between the importance sampling wave function and the exact ground state. The PIMC 
approach of Ref. ^ provides nQ'^'*(PIMC) = 0.069 at temperature T=1.18 K, with large 
statistical errors. As a final comment, we stress that all the above theoretical values of the 
■^He condensate fraction are slightly lower than the latest experimental estimates of Snow et 
a/.,0 r;,Q'^'*(expt) ~ 0.10. However, as the condensate fraction, as well as the kinetic energy, 
is extracted by fitting the Compton scattering profile in neutron scattering experiments at 
large momentum transfers, the resulting ng"^^ can be strongly model dependent. 

We start the analysis of the mixture by studying the Xa-dependence of ^He momentum 
distribution. Fig. 1 shows kn^^\k) / {{271)^ p^) in mixture at X3 = 0.066 (pexpt = 0.358 a^^) 
compared with that of pure ^He (p4 = 0.365 cr^^), both at P=0 . The differences are small 
and can be explained by the slight change in density. In fact, the smaller mass of ^He results 
in a larger zero point motion of ^He compared with ^He, and therefore the total density of 
the mixture decreases when X3 increases. 

Fig. 2 illustrates the same comparison but for the '^He one-body density matrix. The 
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asymptotic value of p^^\r), identified witfi tfie condensate fraction, is reached at r ~ 7 A. 
The value of n'o^ in the mixture is slightly larger than in the pure phase (see also Table I) 
due mainly to the smaller total density of the mixture. The fermionic nature of the ^He 
does not affect Uq^K In fact, one gets the same n^f^ in the boson-boson approximation, which 
consists in treating the ^He component as a bosonic mass-3 one. Furthermore, if ACA is 
assumed, the boson-boson approximation yields a no'''* which is exactly the one of pure ''He 
at the total density of the mixture. 

The Fermi statistics makes the Xa-dependence of n'^^\k) more sizeable. The ^He momen- 
tum distributions at = 0.066 and = 0.020 are compared in Fig. 3. The corresponding 
Fermi momenta are kp — 0.235 and kp — 0.347 A""^, to be compared with kp — 0.79 
A-i for pure ^He at equilibrium density. The Fermi momentum and the discontinuity Zp in- 
crease along xs , whereas the depletion decreases (see Table I) . This behavior is qualitatively 
explained by considering the change of both the total and partial ^He densities. 

p(^) (r) at Xs — 0.066 is compared in Fig. 4 with the free fermionic case {ypir) / p = l{kpr)) 
and with that of pure '^He at the same p^. In this density region it is necessary to reach 
large r- values before p^^\r) begins to oscillate around zero. Despite of the small partial ^He 
density, p^^\r) is very different from those obtained both in the pure (short-dashed line) and 
the free (long-dashed line) cases. While the pure "^He shows a density matrix very similar 
to the free case, the mixture p^^^ (r) has a strong depletion due to the correlations with the 
^He atoms. This behavior translates into a correspondingly large depletion of n^^\k) at the 
origin. The three density matrices have the nodes approximately at the same points, the 
location of the zeros being governed by the zeros of l{kpr). In fact, by taking the lowest 
order term of the expansion of p^^\r) in powers of the statistical correlation l{kpr), as done 
in the Wu-Feenberg expansion for the distribution function, one gets 

p'Ur)=pf(r)'^ (30) 

where Pb\i^) is the '^He density matrix in the underlying boson-boson mixture. Due to the 
small values of x^ in the mixture, Pwp{i^) is almost indistinguishable from the exact p^^^(r). 
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Eq. (pOD explicitly decouples the statistical and dynamical correlations contributions to 
p^^\r) and has also recently proved to describe quite accurately even the pure ^He density 
matrix.i3 In this approximation, n^^\k) is given by 

n^wUk) = r ^'^'""bH k - k' I) . (31) 

[ZTTj-^PS Jo 

Therefore, the discontinuity Zp coincides with the value of the condensate fraction associated 
to n^si^)- The kinetic energy associated to n^)p{k) can be expressed as 

where Tbs/N^ is the kinetic energy associated to (k). In ACA, the density matrices of 
the two components of the underlying boson-boson mixture are the same and equal to the 
density matrix of pure ^He considered at the total density of the mixture. As a consequence, 
the corresponding condensate fractions are also equal and in this model Zp and nl^^ coincide. 

More detailed information on the Xs-dependence of the condensate fraction, the disconti- 
nuity of n^^\k) at the Fermi surface and the kinetic energies of the two components is shown 
in Table I. Ts{x3 = 0.) is the kinetic energy of one ^He impurity in ^He. Recent DMcii 
and PIMci calculations predict a smaller T^ix^ = 0.) value of about 17.5 K. The effect 
of the three-body correlations is similar to that in the ^He pure phase, i.e., they slightly 
decrease the condensate fraction and simultaneously decrease by about half a Kelvin the 
total kinetic energy. The condensate fraction no"*"* shows a small increment with x^. As we 
have mentioned before, this is mainly a consequence of the fact that the total density of the 
mixture slightly decreases when x^ increases. The effect of the Fermi statistics on ti'q'^ is 
almost negligible, the results of nl^^ in the boson-boson approximation being equal to the 
ones reported in Table I. 

n^Q^ is shown in Fig. 5 as a function of the pressure, P, for pure ^He (diamonds) and for a 
Xs = 0.066 mixture (circles). The condensate fraction, in both cases, decreases with pressure 
as a consequence of the corresponding increase of density. The density of pure ^He is larger 
than the one of the mixture at the same pressure and therefore the condensate fraction in 
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the mixture is larger than in ^He. However, as P increases, the differences between the 
densities become smaller and the condensate fractions of both systems get closer. 

The low values of Zp imply a large value of the energy- dependent effective mass at the 
Fermi surface, 

Me = 1- ^^^{P, E) \E=eF,V=PF = (33) 

where E) is the self-energy of the ^He atoms in the mixture. At X3 = 0.04, Me = 12m3 
which is around three times larger than for pure He at the saturation density.Sil This large 
value of the energy- dependent effective mass can be attributed to the correlations with the 
■^He atoms, and implies a small value of the k dependent effective mass in order to reproduce 
the total effective mass that, at those small concentrations, can be taken m^/ms = 2.3,00 
i.e., the value in the impurity case. 

Fig. 6 shows n^'^\k)/p4 and un^^^ / for a 6 % mixture (solid and long-dashed lines re- 
spectively) together with n^^\k)/p4 for pure "^He at the equihbrium density (short-dashed). 
The three momentum distributions are very close above kp, as the large-A; behavior is essen- 
tially dominated by the short-range dynamical correlations. As in the pure phases, the tails 
of the momentum distributions {k > 3.5 A^^) are taken to have an exponential behavior. 
Their contribution at x = 6.6% to the total kinetic energy is ~ 8%. On the other hand, the 
kinetic energy of the free Fermi sea (that would give an upper-bound to the contribution 
to below kp) is 0.58 K. That means that more than 97% of the ^He kinetic energy 

comes from momenta above kp, clearly showing the importance of the correlations between 
^He and ^He atoms. 

It is also of interest to consider the dependence of on the concentration. Fig. 7 

gives in function of the ^He partial density in the mixture along the P = isobar. 

Obviously, the kinetic energy ends up with the kinetic energy of pure ^He (~ 12 K) which 
corresponds to a density value that lies out of the plot. Therefore the kinetic energy of the 
■^He should be in average a decreasing function of the concentration except for the behavior 
at the origin where the term associated with the free Fermi kinetic energy dominates the 
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overall decreasing behavior driven by the decrease of the total density. Actually, the kinetic 
energy in the interval considered here is well parameterized as the sum of the free Fermi gas 
energy plus a linear term describing the decrease of the kinetic energy with the density 

The numerical value of the parameter A may be estimated by calculating the de- 
pendence of the kinetic energy in the underlying boson-boson mixture and it results to be 
A = 27.2 Ka^. 



III. DISCUSSION AND CONCLUSIONS 

The results obtained in this paper for the '^He condensate fraction and the dependence 
of the '^He kinetic energy are in contrast with recent experimental estimates. In fact, Sokol 
et a/.,i0, analyzing deep inelastic neutron scattering measurements carried out for a 9.5% 
mixture at 1.4 K, and for a momentum transfer as high as 23 A~^, estimated a condensate 
fraction nl^^ = 18% and a '^He kinetic energy of approximately 10 K, basically independent 
on the concentration. These results are to be compared with the theoretical predictions 
ti'q^ ~ 10% and ~ 19 K obtained in AC A for a similar mixture. 

It has been arguedi that the main source of discrepancy with a preliminary presentation 
of the present result Jl3 is due to the use of ACA, implying the same type of local environment 
for the different types of atoms in the mixture. Sokol's observation is physically founded on 
the large zero point motion of the ^He atoms which should decrease the local density around 
them to a value similar to the pure ^He. Obviously, the use of optimal correlations should 
clarify this point. However, it must be stressed that the T = DMC calculations of Ref. 
p9| give for the ^He impurity kinetic energy T3 = 17.5 K, i.e. a 1.5 K lower value than the 
ACA prediction estimated by using the pure He DMC kinetic energy (T4 = 14.3 K).0 On 
the other hand, the predicted no"*"* by DMC^ points to an extrapolated value of 11% for a 
6.6 % mixture at the same temperature. A dramatic change of both n'o'^ and T3 at higher 
concentrations would be required in order to reproduce the experimental estimates. 
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In conclusion, we believe that although the use of optimal correlations will certainly 
decrease the kinetic energy of the ^He component and enhance a little the '^He condensate 
fraction, the resulting values will be far from the present experimental analysis. A full 
theoretical calculation of the scattering process including final state interactions and the 
experimental broadening, similar to the ones performed in pure ■*He,@ is necessary in order 
to fully understand the experimental measurements and reliably extract kinetic energies and 
condensate fractions.. 

Summarizing, we have calculated the momentum distributions of ^He-'^He mixtures in 
the framework of the HNC/FHNC equations using variational wave functions with two- and 
three-body correlations. These momentum distributions can be used as input for the analysis 
of the recent performed inelastic neutron scattering experiments. It has been found that, at 
the low concentration where the mixture is stable, the Fermi statistics do not significantly 
modify the value of the ^He condensate fraction. On the other hand, it is crucial to take into 
account the Fermi statistics for the stability of the mixture. The concentration dependence 
of the different quantities studied in the paper can be mainly explained by the decrease in 
the total density of the mixture when the ^He concentration increases. 
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APPENDIX: 

In this Appendix we present the HNC/FHNC equations for the mixture one-body density 
matrices. 
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The sums of the nodal diagrams contributions, N^;^^^^ and Nj;"^, are obtained by solving 
the integral equations 

Ni'-J = E PA E(^i"/^ - Nt'^ - I 9^yf^ - ^y^) > (Al) 

A=3,4 z,y 

and 

Nj^L = P^X9..c + l{hpT^2)ly - ivS I gc.. + II (A2) 

+P3(-V^ I 2(^?e.. + //// - iV£) ) - (^7cc + llv - iVcc)) . 

The notation (A(rjj) | B{rjk)) stands for the convolution product 

(A(r,,) I B{rjk)) = I dr,A{r,^)B{r^k) ■ (A3) 

The summations over z and y (where z,y = d, e, c) always extend to all possible connections 
allowed by the diagrammatic rules of the HNC/FHNC theory.ffl 

Besides the distribution functions giy'^\r) {gdd'^\ gd^'^\ gil'^^ and g^'^^), which have been 
defined elsewhere,00 it is necessary to introduce the auxiliary distribution functions: 

gl^/\r) = f-'^\r)exp[B;^/\r)], (A4) 
gt'\r) = glf \r)Bt''\r) , (A5) 
gS;?ir)=gSf\r)^, (A6) 



where 



Bt'^ (r) = Ar("/) (r ) + (r ) + Cfr/^ (r ) , ( A7) 



and 



i^.(r) = -/(M + z/5£f(r) . (A8) 
The functions El''/\r), E^f/\r) and E^^f (r) give the contributions of the elementary dia- 



grams. 
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The nodal functions N^^'^\r) are solutions of the following integral equations: 

A=3,4 2,2/ 

N^^:? = P3i9^.c - iV^.c + l/iy I 9cc) . (AlO) 
Finally, the functions C^^"'^^ (r) give the contribution of the dressed triplet correlations, 
CL'-fHru) = E PA / rfr3.;('^'^'^)(n2,ri3,r23)E^?i"/Hri3)^7£.''^Hr32) , (All) 

A=3,4 zy 

and 



C^Sf (n2) = P3y rfr3^(3'3'3)(^^^^^^^^^^^)£3,3)(^^^)^(3,3)^^^^) _ (^12) 

The functions N^y'^^^r) and C^y'^^(?") have been defined in Ref. 

The quantities F^"^ and F^"'', entering the expressions of the strength factors no"^ are 
given by 

ri"^ = E PA / rfr(^f/)(r) - 1 - Nir\r) - Eit\r)) 

A=3,4 '' 

+P3 / t/r(^?(^'3)(r) - Nt'\r) - Et^\r)) 

-(1/2) E PA / dr{gi^'\r)-l + 6U':''\r)){Nif\r) + 2Ei-/^ 

A=3,4 

-(l/2)p3 / dr{gif\r) - l)(A^.:'=^)(r) + 24:'3)(r)) 
-(1/2) E PA / rfr(^i"/^(r) + 5,3^?(:'^)(r))cS-^)(r) 

A=3,4 

-(l/2)p3 / drgif\r)Cir\r) + 4") (A13) 



where E^"^ is the sum of the one-point elementary diagrams.i'i0 By setting ps = (p4 = 0), 
expression (2.15) reduces to the pure phases Fj,.!'! 
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TABLES 



TABLE I. ^He condensate fraction, ^He Zp factor and partial kinetic energies in the mixtures 
as a function of the ^He concentration at zero pressure. The first fines are tfie Jastrow values. Tfie 
second lines include the effect of the triplet correlations. 
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FIGURES 

FIG. 1. Momentum distribution of the ^He atoms in the mixture. The continous line cor- 
responds to X3 = 0.066 (p = 0.3582(7^'^) and the dashed Une to pure ^He at saturation density 
(p = 0.365cr~3). gQi-j^ results are at zero pressure. 

FIG. 2. One-body density matrix of the ^He atoms in the mixture. The notation is the same 
as in Fig. 1 

FIG. 3. ^He momentum distributions in the mixture at X3 = 0.066 (solid line) and X3 = 0.02 
(dashed line). The values oi kp are 0.347 and 0.235 A~^, respectively. 

FIG. 4. One-body density matrix of the ^He atoms in a X3 = 0.066 mixture (solid line) com- 
pared with the free Fermi system ( dash-dotted line) and pure ^He (dashed line), both at the same 
partial density p^. 

FIG. 5. Condensate fraction as a function of pressure. The diamonds and circles correspond to 
pure ^He and to a 2:3 = 0.066 mixture, respectively. The lines are guides to the eye. 

FIG. 6. Momentum distributions per particle of pure ^He at equilibrium density (short-dashed), 
and of ^He (long-dashed) and ^He (solid) of a X3 = 0.066 mixture. 

FIG. 7. ^He kinetic energy as a function of at P = 0. The solid line is the fit provided by 
Eq. (H). 
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